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information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

January  2006 _ Conference  Paper  Preprint _ 

4.  TITLE  AND  SUBTITLE 

INVESTIGATION  OF  THE  SELF-PUMPED  TWO-BEAM  COUPLING  IN  A 
PHOTOREFRACTIVE  MATERIAL  USING  BEAM  PROPAGATION 
SIMULATION  (Preprint) 

6.  AUTHOR(S) 

D.R.  Evans,  G.  Cook,  J.  L.  Cams,  and  M.  A.  Saleh,  (Agile  Filters  Project,  Exploratory 
Development) 

P.  P.  Banerjee  (University  of  Dayton) 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Agile  Filters  Project,  Exploratory  Development  University  of  Dayton 

Hardened  Materials  Branch,  Survivability  and  Sensor  Materials  Division  Dayton,  OH 

Materials  and  Manufacturing  Directorate,  Air  Force  Research  Laboratory 
Wright-Patterson  Air  Force  Base,  OH  45433-7750 

_ Air  Force  Materiel  Command,  United  States  Air  Force _ 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 
Materials  and  Manufacturing  Directorate 
Wright-Patterson  Air  Force  Base,  OH  45433-7750 
Air  Force  Materiel  Command 
United  States  Air  Force 


17.  LIMITATION  18.  NUMBER 
OF  ABSTRACT:  OF  PAGES 

SAR  14 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

Dean  R.  Evans,  Ph.D. 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

(937)  255-4588 


16.  SECURITY  CLASSIFICATION  OF: _ 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 


5a.  CONTRACT  NUMBER 

IN-HOUSE _ 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

62102F _ 

5d.  PROJECT  NUMBER 

4348 _ 

5e.  TASK  NUMBER 

RG _ 

5f.  WORK  UNIT  NUMBER 

M03R1000 _ 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-ML-WP-TP-2007-501 

10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/MLPJ 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-ML-WP-TP-2007-501 


l 


Investigation  of  the  self-pumped  two-beam 
coupling  in  a  photorefractive  material 
using  beam  propagation  simulation 


M.  A.  Saleh8,  P.  P.  Banerjeeb, 
G.  Cook3,  Jenny  Cams3,  and  D.  R.  Evans3 

flArr  Force  Research  Laboratory, 
Wright-Patterson  Air  Force  Base,  OH 

University  of  Dayton,  Dayton,  OH 


V 

Outline 

>  Motivation  of  this  work 

>  Beam  propagation  method 

>  Simulation  results 

>  Conclusions 


/<*> 


p(x) 


Ep) 


Ants) 


Jr 

, JL_ L 

/ 

[  j 

Two  plane 
produce  the 
intensity. 


Waves  Interfere  to 
spatially  modulated 


» cm; 


;ei 


The  general 
diffuse  to 
varying  cha 
in  turn  gives 
field  distribgt 
phase  shift' 

The  refractiv 
produced  thi 
optic  effect. 

Can  be  sum; 
simplified  e: 


ed  free  charge  carriers 
ate  the  spatially 
ilge  distribution,  which, 
rise  to  a  electrostatic 
ion  that  is  90-degree 


e  index  variation  is 
rough  the  linear  electro* 


marlzed  by  this 
xpression1 

*Vl{x,y,z) 


J.  J.  Liu  and  P.  P.  Banerjes.  J.  Opt.  Sac.  Am,  B.  111  pp.  16BS-1693,  1994 


>  In  addition  to  the  diffusion  process,  some  anisotropic  crystals  (Le.,  LiNb03) 
exhibit  strong  photovoltaic  effect  which  refers  to  the  tendency  of  the  photo* 
ionization  process  (under  uniform  illumination)  to  eject  the  electron  In  a 
preferred  direction  due  to  the  asymmetric  potential  in  the  anisotropic  crystal. 

A n  =  cpvEs  ccl(x,y,z ) 

>  Resulting  current  flow  can  be  opposed  by  applying  an  external  E-field. 

>  In  the  closed  circuit  condition  the  photocurrent  will  equal  zero  at  -Es ,  the  effects 
of  the  PV  and  ext  E-fieid  cancel  out 

>  Under  open  circuit  condition  it  is  generally  assumed  that  the  bulk  field  cancels 

out  the  photovoltaic  field  and  PV  effect  does  not  assist  TBC. _ 


V  Contra-directional  Self-pumped  JKt 
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Periodic  index  modulation  from  the  phoioref r active  effect  can  be  used  to 
couple  energy  between  two  interfering  beams*  One  novel  variation  of  the 
contra-directional  beam  coupling  is  the  self-pumped  beam-coupling  where 
the  signal  beam  is  generated  by  the  Fresnel  reflection  of  the  pump. 


>  Strong  dependence  of  the  TBC  efficiency  on  the  spot  size  Is  observed, 
attributed  to  diffraction  and  the  resulting  change  in  the  spatial  overlaps. 

>  The  goal  of  our  simulation: 

*  Is  to  simulate  arbitrary  shaped  beam-coupling  using  exact  diffraction. 

*  is  to  analyze  the  contribution  from  photovoltaic  effect  and  diffusion  in 
two-beam  coupling. 

>  To  the  best  of  our  knowledge,  self-pumped  contra-directional  TBC  has  not 
been  done  with  BPM. 
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Beam  Propagation 
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EJy)  as  sum  of  discrete  plane  waves 
of  different  amplitudes 


Spatial  amplitude  distribution  is  expressed  as  superposition  of  plane  waves. 
Each  plane  wave  accumulates  different  phase  going  from  one  z-plane  to 
anther  z-plane.  After  a  distance  Az  accumulated  phase  is  added  to  each  plane 
wave  and  the  spatial  amplitude  distribution  can  be  recalculated  again. 
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The  optical  path  is  broken  into  a  sequence  of  finite  steps. 


Simulated  Interference 
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Steps  of  simulation 


>  We  simulate  the  interference  of  the  beams  by  adding  the  forward  and 
backward  propagating  beams  at  every  incremental  z  value.  This  produces 
the  modulated  /(z)  inside  the  crystal. 

We  use  this  modulated  Intensity  to  calculate  the  space  charge  field. 

Change  in  index  then  can  be  calculated  from  the  electrooptic  effect. 

>  We  modify  the  index  inside  the  crystal  by  using  the  simplified  expression1 
(omitting  the  details  of  space  charge  field  generation): 

n  =  n0+An,  where  A n=cpvEsc<x  I(x,y,z)  (PV) 

or,  A n  =  cdfEsc^Vl(x,y,z)  (diffusion) 

>  For  this  work  only  a  X  Jong  crystal  was  simulated.  We  chose  the  two 
constants  such  that  there  was  observable  beam  coupling  in  this  thin  crystal 
and  the  An  is  same  for  both  diffusion  and  PV  effect. 

J.  J.  Liu  and  P  P.  Banerje*.  J.  Opt  SOC  Am  B.  11,  pp  168S-1693,  1994V 


Condition  for  Beam  coupling 


>  This  total  phase  difference  comes  from  the  phase  difference  between 
interference  pattern  and  grating  pattern  and  from  phase  change  due  to 
the  reflections  from  the  grating. 

>  A  pump  (signal)  power  is  defined  as  the  change  in  pump  (signal) 
transmission  between  linear  case  and  after  grating  formation. 

>  A  pump  power  and  A  signal  power  tells  if  destructive  or  constructive 
interference  has  taken  place. 


Steps  of  convergence 


>  For  a  given  intensity  distribution  J(x,y,z)  inside  the  crystal,  we 
calculate  n (x,y,z) . 

>  We  obtain  a  new  intensity  distribution  /{x,y, z)  through  another 
round  of  beam  propagation,  where  p unrip  and  signal  scatters  are 
also  accounted  for,  which  gives  us  a  new  n(x,y,z). 

>  Error  is  defined  as  the  average  percent  difference  between 
successive  interference  amplitude  arrays. 

>  We  continue  the  process  until  we  obtain  a  converged  steady  state 
solution  for  both  /(x^)  and  n(x#p)  using  iterative  shooting 
method. 

>  For  the  converged  /(xsy, z)  and  n(xj,z)  pump  and  signal  scatters  are 
added  to  signal  and  pump  respectively  and  the  power  changes  due 
to  the  interference  are  calculated. 


Considering  diffusion  alone 


The  index  and  interference  are  90-degree  out  of  phase  (as  expected).  An  =  0.022 
{for  Cm  =  0.27).  AP  =  0,58  %  of  pump  power.  0^=  >0.27  (equivalent  to  rotating 
optic  axis  180° )  makes  AP  =  -0.58  %  of  pump  power. 


Considering  both  PV  and  Diffusion 


Phase  difference  >  0  and  <  90  deg.  An  =  0.031 ,  larger  than 
considering  them  alone.  AP  changes  with  PV  contribution. 


Photovoltaic  Contribution 
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Observed  photovoltaic  contribution 
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G.  Cook,  J.  P.  Duignan,  D.  G. 
Jones,  “Photovoltaic  contribution 
to  counter-propagating  two-beam 
coupling  in  photorefractive  lithium 
niobate"  Optics  communications, 
voi.  192,  pp  393-398,  2001. 


Applied  voltage  (kV) 
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_ E0  +  iED  +  Epy  _ 

1  +  EDtEQ -i(E0/EQ  +(Na/NdXEpv /Eq)) 


Conclusions  and  Future  Work 


>  At  this  point  we  are  capable  of  simulating  the  self  pumped 
contra-directional  two-beam  coupling  in  a  photorefractive 
medium  with  arbitrary  shaped  beams. 

>  Our  simulation  shows  positive  role  of  the  photovoltaic  effect  in 
self-pumped  contra-directional  TBC,  in  agreement  with 
experimental  observation. 

>  A  more  complete  simulation  that  includes  dark  conductivity  and 
un-simplified  space  charge  field  is  in  progress. 


